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Abstract 

We compute the exact number, L{n), of solutions to the Langford 
pairings problem for any positive iirteger n < 29 and the exact num¬ 
ber of solutions to the Nickerson variant of the problem, N{n), for any 
positive integer n < 26. These numbers correspond to the sequences 
r014552, r059106 in Sloane’s Online Encyclopedia of Integer Sequences. 

The exact value of these numbers were known for any positive integer 
n < 27 for the A014552 sequence and for any positive integer n < 24 
for the A059106 sequence. First we report that the number of Langford 
pairin g s for n = 27 is L(27) = 111,683,611,098,764,903,232, and for 
n = 28 it is L(28) = 1,607,383,260,609,382,393,152. Next we report 
that the number of solutions for the Nickerson variant of Langford pair¬ 
ings for n = 24 is iV(24) = 102, 388, 058, 845, 620, 672 and for n = 25 it is 
A(25) = 1, 317, 281, 759, 888, 482, 688. 

1 Historical Background 

The Langford problem is named after C. D. Langford, a Scottish mathematician, 
who devised the problem in 1958 after having observed his son playing with 
colored blocks. He noticed that his son had arranged a set of six block, two 
red, two blue, and two green, in such away that the pair of red blocks were 
separated by a single block, the blue pair by two blocks, and the green pair by 
three blocks. He further noticed that he could add a pair of yellow blocks to the 
arrangement in away that would preserve the distances of the previous blocks 
while having the yellow pair separated by four blocks. 



He captured this idea using numbers and asked the question. Given a sequence 
of 2n numbers {1,1, 2, 2,..., n, n} find a permutation in which the two copies 
of each number k are k units apart. For instance, when n = 4 the permutation 
2, 3,4, 2,1, 3,1,4 would be a solution. 

At about the same time while working on Steiner triple systems, Norwegian 
mathematician T. H. Skolem, in 1957, proposed a similar problem, he asked if 
it was possible to distribute the number {1,2,..., 2n} in n pairs (or, b^) such 
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that br — ttr = r for r = 1,2,..., n ?”. For instance when n = 4, the pairings 
(1,4), (2, 6), (3, 5), (7, 8) is one solution. 

In 1966 Nickerson, unaware of Skolem’s problem, proposed a varient of Lang¬ 
ford’s problem where the pair of numbers k are separated by exactly k—1 other 
numbers. For instance when n = 4 the sequence 3,4, 2, 3,2,4,1,1 is a solution. 
Notice that if we consider the sequence as being placed in an array with 2n 
positions and index from 1 to 2n then this sequence is the same as the Skolem 
pairing. 
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The first natural question to asked is when do such pairings exist, that what 
are the necessary and sufficient conditions for the existance such sequences. In 
1959 Daveis answers this question for the Langford formulation, and showed 
the sequence exists whenever n = 0,3 mod 4, following the observation that 
that the pairs of number occupy the positions 1 to 2n in some order, the first 
occurrence of k is at position Ok and the second at -I- fc -f 1. Furthermore he 
gave a method of construction such a sequence. 


1—1 

1 

II 

4k-4, ..., 2k, 4k-2, 2k-3, ..., 1, 4k-l, 

1, ..., 2k-3, 2k, ..., 4k-4, 2k-l, 

4k-3, ..., 2k+l, 4k-2, 2k-2, ..., 2, 

2k-l, 4k-l, 2, ..., 2k-2, 

2k-fl, ..., 4k-3 
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4k-4, ..., 2k, 4k-2, 2k-3, ..., 1, 4k-l, 

1, ..., 2k-3, 2k, ..., 4k-4, 4k, 

4k-3, ..., 2k~fl, 4k-2, 

2k-2, ..., 2, 2k-l, 4k-l, 2, ..., 2k-2, 

2k-(-l, ..., 4k-3, 2k-l, 4k 
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Figure 1: Davies construction of a solution to the Langford pairing 


In the case of the Skolem problem the existance question was answered by 
Skolem himself in 1957. He showed that the pairings exists whenever n = 0,1 
mod 4, following a similar observation that pairs of number occupy the positions 
1 to 2n in some order, the first occurrence of k is at position Cfc and the second 
at Ofc -l- k. And a method of construction such a sequence was presented as well. 


n = ‘ik 


(4fc -I- r, 8A: — r) for r = 0,1,..., 2fc — 1, 
(2fc -f l,6fc), (2fc,4fc - 1), 

(r, 4fc — 1 — r) for r = 1, 2,..., fc — 1, 

(k, fc -I- 1), 

(fc + 2 -f r, 3fc — 1 — r) for 0,1,..., fe — 3. 



n = 4fc -f 1 


(4fc -I- 2 -I- r, 8fc + 2 — r) for r = 0,1,.. ., 2fe — 1, 
(2fc -f 1, 6fc + 2), {2k -f 2, 4fc + 1), 

(r, 4fc + 1 — r) for r = 1,2,... ,k, 

{k l,k 2), 

{k + 2 + r,3k + 1 — r) for r = 1, 2,..., fc — 2. 



Figure 2: Skolem’s construction of a solution to the Skolem pairing 
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A014552 I A059106 | A176127 | A004075 | 
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26 

0 

52 

0 

8 

150 

252 

300 

504 

9 

0 

1 , 328 

0 

2 , 656 

10 

0 

0 

0 

0 

11 

17 . 792 

0 

35.584 

0 

12 

108,144 

227 , 968 

216,288 

455,936 

13 

0 

1 , 520,280 

0 

3 , 040,560 

14 

0 

0 

0 

0 

15 

39 . 809,640 

0 

79 . 619,280 

0 

16 

326 , 721,800 

700 . 078,384 

653 . 443,600 

1 , 400 , 156,768 

17 

0 

6 , 124 , 491 , 248 

0 

12 , 248 , 982,496 

18 

0 

0 

0 

0 

19 

256 , 814 , 891,280 

0 

513 , 629 , 782,560 

0 

20 

2 , 636 , 337 , 861,200 

5 , 717 , 789 , 399,488 

5 , 272 , 675 , 722,400 

11 , 435 , 578 , 798,976 

21 

0 

61 , 782 , 464 , 083 , 584 

0 

123 , 564 , 928 , 167 , 168 

22 

0 

0 

0 

0 

23 

3 , 799 , 455 , 942 , 515,488 

0 

7 , 598 , 911 , 885 , 030,976 

0 

24 

46 . 845 , 158 , 056 , 515,936 

102 , 388 , 058 , 845 , 620 , 672 

93 , 690 , 316 . 113 , 031.872 

204 , 776 , 117 , 691 , 241 , 344 

25 

0 

1 , 317 , 281 , 759 , 888 , 482 , 688 

0 

2 , 634 , 563 , 519 , 776 , 965 , 376 

26 

0 

0 

0 

0 

27 

111 . 683 . 611 , 098 , 764 , 903 , 232 * 

0 

223 . 367 , 222 , 197 , 529 . 806 , 464 * 

0 

28 

1 , 607 , 383 , 260 , 609 , 382 , 393,152 

7 7 ? 

3 , 214 , 766 , 521 , 218 , 764 , 786,304 

7 77 


Table 1: Langford number 4.014552, Nickerson’s variant 4059106, 4014552 = 
2 X 4014552, and Skolem number 4059106 = 2 x 4059106 for n > 1, * n = 27 
was previously reported to be 111, 683,606, 778,027,803,456 


In this note we are concered with the problem of counting the number solution. 
In what follows we present a breif overview of algorithms for counting these 
sequences, implementation and new results we found. 


2 Counting Algorithms 

There are in general three approaches to counting the number of pairings. First 
the generate and count methods, second the algebraic method, and third a 
method based on counting technique of inclusion exclusion. All the procedures 
below will be given interms of counting the number of Langford pairing, counting 
all solutions where reflected solutions are counted as distinct. 

2.1 Miller’s Generate and Connt 

One method of counting is to generate all possible permutations of the 2n num¬ 
bers in some order and count as each valid solution is generated. However this 
is quit inefficient since there are permutations. A more efficient procedure 
is to use a backtracking algorithm for generating the solutions in a systematic 
way and counting them. One such procedure is given by Miller on his website 
[1], and reproduced below. The alorithm proceeds by placing the pairs in de¬ 
creasing order, starting with leftmost available position where the pair can fit 
into. Once the pair is place it will try to palce the next smaller pair if it can 
not be placed then the previously placed pair needs to be moved to the next 
available valid position. The algorithm stops when all possible positions for the 
largest pair have been explored. 

• Let n be the number of pairs to be placed 
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• Let A be an array with 2 n positions holding the current (partial) arrang- 
ment 

• Let k keep track the pair of numbers that are to be placed into A 

• Start by inserting the highest pair (k = n) into the first position. The 
pair will fit because the arrangement is empty. 

• Next, decrement k and find the left most position p where both A[p] and 
A[p + fc + 1] are empty 

• If such a position p is found, then let A[p] = A[p + fc + 1] = fc. 

— If the I’s were just placed, a solution is reached, count it, and re¬ 
move the I’s, 2’s and 3’s from A^ and continue scanning for the next 
available placment of fc = 3. 

• If the fc*^ pair can not be placed, forget about placing them, and increment 
the value of fc. remove fc’s from A, and find a new placement for the fc*^ 
pair starting from the next position to the right. 

~ If fc is incremented to beyond n, then the all the solutions have been 
found and algorithm terminates 


2.2 Godfrey’s Algebraic Method 

An algebraic method for computing the number Langford pairings was proposed 
by Godfrey. There is no ofhcial paper on the algorithm but the description of this 
method appears in [9]. The problem is modeled by a polynomial where each term 
represents the label and its position, and the number of pairings is then the co¬ 
efficient of the term xiX 2 ■ ■ ■ Xn- Consider L(3), and X = {x:i^X 2 ,xz, 2 : 4 , X 5 ,xe), 
then 


F(3, X) = {xixs + X 2 X 4 + X 3 X 5 + X 4 Xe){xiX 4 -I- a;2a;5 -I- a;3a;6)(a;ia;5 -I- a;2a;6) 

each of the factors represents the possible ways in which a label and position 
can appear in the solution. For instance 3 can appear either in the first and fifth 
positions or in the second and sixth. When the polynomial is expanded then for 
example a term like {x 2 X 4 ){x 3 Xq){xiX 3 ) correspond to a sequence 3,1,2,1,3, 2, 
and the coefficient to the term a;ia; 22 ; 3 a; 4 a; 5 a ;6 will be the number of possible 
pairings (twice that since symmetric solutions are considered the same). 


F(3, X) = xJa:3X4i:5 + x^X4^x^xq + xiX2X^x^ + 


Xg + 3 
X ? + 3 


^6 + 
Xfl + X 


1x214x5 - 


C2X3X4Xg + x^xgxgxg + 3 

■3X5X6 + x|x 3 X 4 X 6 + x|x 


lX2X3X4xe + xiX2X3X5xe + xiX2X4X5xe 


2,22 , 
X2X4X5 + X1X4X5X6 + 2 

XgXg + X1X3X5X6 + XIX 

+ X1X3X4X5X6 + 2 xix - 


However computing the coefficient by expanding the polynomial is just as time 
consuming as the generate and count methods since there are ^‘^^^ 2)1 ■ over 
come this the variables xi,, X 2 n are allowed to take on the values 1 and — 1, 
and the resulting values of {xi x ... x X 2 n) x F{X, n) are summed, it is easy to 
see that the result is 2^"L(n) (twice that incase symmetric solutions are counted 
as the same) 
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( 2n \ n 2n—i—l 

) n X! XkXk+i+i = 2 ^”L(n) 

j=l / i=l fc=l 

Each of the terms in F other than X 1 X 2 ■ ■ ■X 2 n is missing at least one of the 
variables (ife), then these terms when multiplied by Xk result in zero after sum¬ 
ming over Xk = +1 and —1. The calculation costs 2^” evaluations of F{n,X), 
each evaluation takes O(n^) multiplications and additions. 

In practice the cost of multiplications can be reduced to 0{n) by going through 
the { — 1,1}^” strings in Gray code order, there by replacing the multiplication 
by an addition and subtraction. 


2.3 Larsen’s Inclusion Exclusion Method 


In 2009 Larsen [8] proposed a new algorithm based on the principle of inclusion 
exclusion to count the number of pairings. Inclusion exclusion is a counting 
technique for obtaining the number of elements in the union of sets. In its 
general form the principle of inclusion exclusion for finite sets Ai,..., An states 
that the number of elements in the union of the sets is the sum of cardinalities 
of set intersections, subtracting when the number of sets is even and adding 
when its odd. 








It is common to generalize this form to count the number of elements that are 
not in any of the subsets of Ai. This complementary form is obtained by letting 
5 be a finite universal set containing all of the Ai and letting Ai denote the 
complement of Ai in S', then 


2=1 


n 


n 


s-U^, 

2=1 


or 



i: (-1)''' n-^* 

/C{l,...,n} iG/ 


where ^0 = S since intersecting with no sets will result in all the elements. 
This is connected to the number of pairings by viewing each set Ai as the set 
of (invalid) sequences that avoid the position i. Then Ai n ... n A„ is the set of 
sequences avoid no positions, and a valid solution to the problem. Then L(n) 
can ne expressed in the follwoing form 


L{n) = E 

XC{l,...,2n} 


where X denotes the positions to be avoided and a{X) is the function that 
counts the set of solutions that avoid the positions in X. Let X be a binary 
sequence (xi,... ,X 2 n) G {0,1} where every 0 in the sequence is a disallowed 
poisition. Then 

n 2n—fc—1 

a{X) = ^ ' ^j^i+k+i 

fe=i j=i 
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To compute L{n) need to sum over all binary strings of length 2n. 


n 2n—k—l 

Hn)= ((-1)^*='''') n J2 XjXj+k+1 

{xi,...,x2n)e{o,i} fc=i j=i 

There will be a total 2^" computations of a(X), where each computation costs 
O(n^) time. The total running time then is 0(n^2^”). 

As before the evaluation cost of a{X) can be reduced by going through the 2^” 
binary strings in Gray code order. The advantage using Larsen’s method over 
Godfrey’s is that there is no longer an addition factor 2^”. However Larsen’s 
experiments show that Godfrey’s procedure is faster in practice. 


3 Results 

We implement Godfrey’s procedure for computing the number solutions to the 
Langford and the Nickerson’s variant on NVIDIA’s CUBA parallel comput¬ 
ing platform. The calcualtion is done modulo primes then the actual value 
is obtained by Chinese remainder theorem. We checked or implementation by 
computing the previously known values. New results are found for the original 
problem when n = 27, n = 28. We note that recently the 27*^ langford number 
was reported, and does not match with our results. It is worth noting that our 
results match in the order of magituide and the seven most significant digits. 
We also found new results for the variant when n = 24, n = 25, and partially 
computed the value for n = 28. These correspond to sequences A014552 and 
A059106 in Sloane’s Online Encyclopedia of Integer Sequences. Additionally 
the corresponding values for the sequence of number A176127 and A004075 in 
the oeis are just twice that of A014552 and A059106 respectively, since in the 
orignal problem refelected solutions were considered the same. The table below 
shows the sequences along with our results in red. 
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